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1 Introduction 



Event topologies with a single jet with large transverse energy and large missing transverse 
momentum, referred to as monojets in the following, are important final states for searches 
for new phenomena beyond the Standard Model (SM) at a hadron collider. The primary 
SM process that results in a true monojet final state is Z-boson production in association 
with a jet, where the Z boson decays to two neutrinos. A further important reducible 
contribution to this final state consists of events that include a W boson and a jet, where 
the charged lepton from the TV-boson decay is not reconstructed. 

Phenomenological scenarios beyond the Standard Model (BSM) that result in a mono- 
jet final state include supersymmetry [1-11] and large extra dimensions (LED) [12]. A 
model-independent treatment of the production of dark matter (DM) particles at the Large 
Hadron Collider (LHC) has been proposed recently, where DM particles are pair-produced 
in association with a jet [13-15]. In the following, a search for an excess of monojet events 
over SM expectations is performed. The results are interpreted in a framework of LED and 
DM particle pair production. They are based on a dataset of 4.7 fb _1 of proton-proton 
(pp) collisions at \fs = 7 TeV from the LHC and supersede those presented in the 2010 
ATLAS monojet analysis that used 35 pb _1 of data [16]. Other monojet searches were 
performed in Run I and Run II at the Tevatron [17-19] and also by CMS with the 2010 [20] 
and 2011 [21] LHC datasets. None of these found evidence of new phenomena beyond the 
Standard Model. 

Models of large extra spatial dimensions have been proposed to remove the hierarchy 
problem [22-25] by addressing the weakness of gravity relative to all other forces. One 
popular model of LED that is often used to interpret the results of monojet searches at 
particle colliders is that of Arkani-Hamed, Dimopoulos, Dvali (ADD) [12]. In this model, 
gravity propagates in the (4 + n)-dimensional bulk of space-time, while the SM fields 
are confined to four dimensions. The large observed difference between the characteristic 
mass scale of gravity (the Planck mass) and the electroweak scale (as characterised by the 
W^-boson mass) is the result of the four-dimensional interpretation of the Planck scale, 
Mpi = 1.2 x 10 19 GeV, which is related to the fundamental (4 + n)-dimensional Planck 
scale (M D ) by M Pl 2 = 8vr M D 2+n R n , where n and R are the number and size of the 
extra dimensions, respectively. An appropriate choice of R for a given n results in a 
value of Md close to the electroweak scale. The extra spatial dimensions are compactified, 
resulting in a Kaluza-Klein tower of massive graviton modes. At hadron colliders, these 
graviton modes can be produced in association with a jet. The production processes include 
qg — > qG, gg — > gG, and qq — > gG, where G stands for the tower of gravitons, q for a quark, 
and g for a gluon. As gravitons do not interact with the detector, these processes give rise 
to a monojet signature [26]. 

Particle dark matter is a well-established paradigm to explain a range of astrophysical 
measurements (see for example ref. [27] for a recent review). Since none of the known 
SM particles are adequate DM candidates, the existence of a new particle is hypothesised, 
with properties suitable to explain the astrophysical measurements. One class of parti- 
cle candidates of interest for searches at the LHC consists of weakly interacting massive 
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particles (WIMPs) [28]. These are expected to couple to SM particles through a generic 
weak interaction, which could be the known weak interaction of the SM or a new type of 
interaction. Such a new particle is a cold dark matter candidate, which can be produced at 
the LHC. It results in the correct relic density values for non-relativistic matter in the early 
universe [29] , as measured by the WMAP satellite [30] , if its mass lies in the range between 
a few GeV and a TeV and if it has electroweak-scale interaction cross sections. The fact 
that a new particle with such properties can be a thermal relic of the early universe in ac- 
cordance with the WMAP measurements is often referred to as the WIMP miracle. Many 
new particle physics models designed to solve the hierarchy problem also predict WIMPs. 

Because WIMPs do not interact with the detector material, their production leads to 
signatures with missing transverse momentum (p™ lss ), the magnitude of which is called 
E™ 1SS . Searches involving E™ 1SS at the LHC are therefore canonical WIMP searches, al- 
though the LHC experiments cannot establish whether a WIMP candidate is stable on cos- 
mological time scales and hence a DM candidate. In some super symmetric models, WIMPs 
are expected to be dominantly produced in cascade decays of heavier unstable supersym- 
metric particles along with high-pi SM particles. In a more model-independent approach, 
WIMP pair production at colliders is proposed to yield detectable E™ 1SS if the WIMP pair 
is tagged by a jet or photon from initial- or final-state radiation (ISR/FSR) [13, 31]. Even 
though this approach does not rely on a specific BSM scenario, it does have assumptions: 
WIMPs are pair-produced at the LHC and all new particles mediating the interaction 
between WIMPs and the SM are too heavy to be produced directly; they can thus be 
integrated out in an effective field theory approach. The resulting interaction is hence a 
contact interaction between the dark sector and the SM. It is worth noting that the DM 
particles are not explicitly assumed to interact via the weak force. They may also couple 
to the SM via a new force. Throughout this work, the terms WIMP and DM particle 
(candidate) are synonymous. 



Name 


Initial state 


Type 


Operator 


Dl 


qq 


scalar 


m„ - - 

jfsxxqq 


D5 


qq 


vector 


Tpxi^xqi^q 


D8 


qq 


axial-vector 




D9 


qq 


tensor 


-wixo^xqo^q 


Dll 


99 


scalar 


j^XXa s (G a ^) 2 



Table 1. Effective interactions coupling Dirac fcrmion WIMPs to Standard Model quarks or gluons, 
following the formalism of ref. [32] . The tensor operator D9 describes a magnetic-moment coupling. 
The factor of the strong coupling constant a s in the definition of Dll accounts for this operator 
being induced at one-loop level. G Mt , is the colour field-strength tensor. 

It is assumed here that the DM particle is a Dirac fermion x, where the only difference 
for Majorana fermions would be that certain interaction types are not allowed and that 
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the cross section for each operator is larger by a factor of four. Five interactions are 
considered (table 1), namely Dl, D5, D8, D9, Dll, following the naming scheme of ref. [32]. 
Dl, D5, D8, and D9 describe different bilinear quark couplings to WIMPs, qq — > XX-, an d 
Dll describes the process gg — > XX- The 14 operators for Dirac fermions in ref. [32] fall 
into four categories with characteristic E!f lss spectral shapes. Dl, D5, D9, and Dll are 
a representative set of operators for these four categories, while D8 falls into the same 
category as D5 but is listed explicitly in table 1 because it is often used to convert LHC 
limits into limits on DM pair production. In the operator definitions in table 1, M* is the 
suppression scale of the heavy mediator particles that are integrated out. The use of a 
contact interaction to produce WIMP pairs via heavy mediators is considered conservative 
because it rarely overestimates cross sections when applied to a specific BSM scenario. 
Cases where this approach is indeed optimistic are studied in refs. [15, 33]. The effective 
theory provides a useful framework for comparing LHC results to direct or indirect dark 
matter searches. Within this framework, interactions of SM and DM particles are described 
by only two parameters, the suppression scale M* and the DM particle mass m x . 

2 Data and simulated samples 

The ATLAS detector [34, 35] at the LHC covers nearly the entire solid angle around the col- 
lision point. It consists of an inner tracking detector surrounded by a thin superconducting 
solenoid, electromagnetic and hadronic calorimeters, and an external muon spectrometer 
incorporating large superconducting toroidal magnets. A three-level trigger system is used 
to select interesting events for recording and subsequent offline analysis. Only data for 
which all subsystems described above were operational are used. Applying these require- 
ments to pp collision data, taken at a centre-of-mass energy of -^/s = 7 TeV with stable 
beam conditions during the 2011 LHC run, results in a data sample with a time-integrated 
luminosity of 4.7 fb _1 , determined with an uncertainty of 3.9% [36, 37]. 

Monte Carlo (MC) simulations are used both as part of the background estimation 
and to model signal processes. Processes that dominate the background are Z- or W- 
boson production in association with jets, which are simulated with ALPGEN [38] using 
the parton distribution function (PDF) set CTEQ6L1 [39]. The W -)■ iv plus jets and 
Z — > vv plus jets samples are simulated with up to six additional partons at leading order, 
while the process Z/j* — > £ + £~ plus jets is simulated with up to five additional partons 
at leading order. Additional jets are generated via parton showering, which, together 
with fragmentation and hadronisation, is performed by HERWIG [40, 41]. The MLM [42] 
prescription is used for matching the matrix-element calculations to the parton shower 
evolution. JIMMY [43] is used to simulate the underlying event. Additional Z/W plus jets 
samples generated with SHERPA [44] are used to estimate the uncertainties related to the 
event generator. Single top quark and pair production are simulated with MCONLO [45], 
fixing the top-quark mass to 172.5 GeV, and using the next-to-leading-order (NLO) PDF 
set CTEQ6.6 [46]. Parton showering and hadronisation are performed with HERWIG, and 
JIMMY is again used for the underlying event. Di-boson (WW, WZ, ZZ) samples are 
generated with SHERPA. Backgrounds from QCD multijet production are estimated from 
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data (see section 5.2 ). PYTHIA [47] simulations of this process, normalised to data, are 
used in figures for illustrative purposes only. 

For graviton production in the ADD model, a low-energy effective field theory [26] with 
energy scale Md is used to calculate the signal cross section considering the contribution 
of different graviton mass modes. Signal samples corresponding to a number of extra 
dimensions varying between two and six are considered, with the renormalisation and 



mode produced in the event and pt denotes the transverse momentum of the recoiling 
parton. The samples are produced with an ADD implementation as a user model of PYTHIA, 
which is also used for parton showering and hadronisation. MRST2008 LD** [48] PDF sets 
are used for the event simulation. The event yields for CTEQ6 . 6 PDFs are obtained by re- 
weighting these samples, and are used to estimate cross sections, as well as PDF systematic 
uncertainties. ADD cross sections are calculated at both leading order (LO) and NLO. The 
NLO calculations take into account QCD corrections to graviton production and have been 
produced for the kinematic regions explored here following ref. [49]. 

The effective field theory of WIMP pair production is implemented in MADGRAPH5 [50] 
(version 1.3.33), taken from ref. [32]. WIMP pair production plus one and two additional 
partons from ISR/FSR is simulated requiring at least one parton with a minimum trans- 
verse momentum of 80 GeV. Only initial states of gluons and the four lightest quarks are 
considered, assuming equal coupling strengths for all quark flavours to the WIMPs. The 
mass of charm quarks is most relevant for the cross sections of the operator Dl (see table 1) 
and it is set to 1.42 GeV. The generated events are interfaced to PYTHIA for parton show- 
ering and hadronisation. The MLM prescription is used for matching the matrix-element 
calculations of MADGRAPH5 to the parton shower evolution of PYTHIA. The CTEQ6L1 PDF 
set is used for the event simulation. The MADGRAPH5 default choice fo r the reno rmalisation 
and factorisation scales is used. The scales are set to the sum of \J m 2 + p\ for all pro- 
duced particles. Events with WIMP masses between 10 and 1300 GeV are simulated for 
four different effective operators (Dl, D5, D9, Dll). In all cases, WIMPs are taken to be 
Dirac fermions, and the pair-production cross section is calculated at LO. 

The background MC samples use a detector simulation [51] based on GEANT4 [52] and 
are reconstructed with the same algorithms as the data. The signal MC samples employ a 
mix of the detailed GEANT4 detector simulation and a simulation relying on parametrisations 
of calorimetric signals to shorten the CPU time required (ATLFAST-II [51]). Individual sig- 
nal MC samples have been validated against the more detailed detector simulation relying 
fully on GEANT4. Effects of event pile-up — multiple pp interactions occurring in the same or 
neighbouring crossing of two proton bunches — are included in the simulation. MC events 
are re-weighted to reproduce the distribution of the number of collisions per bunch crossing 
observed in the data. 

3 Analysis strategy and physics object reconstruction 

This analysis is based on reconstructed jets, electrons, muons, and E™ ss . The definitions 
of electron and muon candidates and of E™ 1SS are different when selecting a data sample 



factorisation scales set per event 




where Mq is the mass of the graviton 
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that is used to search for a BSM excess than when selecting a data sample that is used 
to determine a background contribution in the BSM search sample. Such a BSM search 
data sample is referred to as a signal region (SR) in the following and may be sub-divided 
further into a number of potentially overlapping signal regions corresponding to different 
kinematic selection criteria based on Elf lss and p^* 1 , the transverse momentum of the 
highest-pT (leading) jet. A data sample used for the background determination is called a 
control region (CR). The control regions are defined by the physics process used to calculate 
the background. They are further sub-divided into the same kinematic categories as the 
signal regions. Additional selection requirements are imposed in the control regions to 
define data samples enriched in events from certain processes such as W — > /iz^+jets. Such 
a data sample is referred to as the W — > [iv control region. The search for a BSM excess in 
monojet final states is performed by estimating the contribution of background processes 
to the signal regions, either from data control regions or from simulated samples, and then 
searching for an excess of events in the signal regions beyond the estimated backgrounds. 

Jet candidates are reconstructed using the anti-fct clustering algorithm [53] with a 
radius parameter of 0.4. The inputs to this algorithm are clusters of energy deposits in 
calorimeter cells seeded by those with energies significantly above the measured noise [54]. 
Jet momenta are calculated by performing a four-vector sum over these cell energy clusters, 
treating each cluster as an (E, p) four-vector with zero mass. The direction of p is given by 
the line joining the nominal interaction point with the calorimeter cluster. The resulting jet 
energies are corrected to the hadronic scale using px and pseudorapidity (rj) dependent 1 
calibration factors based on MC simulations and validated by extensive test beam and 
collision data studies [55]. 

Electron candidates are required to have pr > 20 GeV and \rj\ < 2.47, in order to be 
within the acceptance of the tracking system. In the monojet signal regions, where events 
with electrons are vetoed, relatively loose criteria are used to define an electron candidate 
(SR electron). A looser electron definition leads to a more stringent veto. SR electrons are 
required to pass the medium electron shower shape and track selection criteria described 
in ref. [56]. No spatial isolation is required. In data control regions used for background 
estimates, events with electrons are explicitly selected. More stringent electron selection 
criteria (CR electron) are used when exactly one electron is selected in order to better 
suppress contamination in these control regions. A CR electron is required to pass the 
tight [56] electron shower shape and track selection criteria mW ^ ei; control regions. In 
addition, the following isolation criterion is imposed for CR electrons to suppress events 
where a jet is mis-identified as an electron: the scalar sum of the transverse momentum 
of tracks with AR = \J (Ac/)) 2 + (An) 2 < 0.2 around the electron candidate, excluding the 
electron itself, has to be less than 10% of its transverse energy (E^). In control regions 
where exactly two electrons are required, the looser SR electron definition without isolation 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) 
in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the 
centre of the LHC ring, and the y-axis points upward. Polar coordinates (r, <j>) are used in the transverse 
(x,2/)-plane, cj> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of 
the polar angle as rj = — lntan(#/2). 
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requirements is used. 

A muon candidate used in the definition of the signal regions (SR muon) is recon- 
structed either by associating a stand-alone muon spectrometer track with an inner detec- 
tor track, or from an inner detector track that is confirmed by a directional segment in 
the muon spectrometer [57]. SR muons, which are used in the veto for signal regions, are 
required to have p-r > 7 GeV and \rj\ < 2.5. They are also required to be isolated: the 
scalar p^ sum of tracks within AR = 0.2 around the muon track, excluding the muon itself, 
must be less than 1.8 GeV. As for electrons, the muon selection criteria in control regions 
are more stringent. A CR muon candidate must have a stand-alone muon spectrometer 
track associated with an inner detector track. Those SR muons that have only an inner de- 
tector track tagged by a segment in the muon spectrometer do not satisfy the CR selection 
criteria. Furthermore, CR muons satisfy > 20 GeV and \r]\ < 2.4, and have an impact 
parameter along z with respect to the reconstructed primary vertex of \zq\ < 10 mm to 
reject cosmic-ray muons. The CR muons are also required to be isolated: the scalar p^ 
sum of tracks within AR = 0.2 around the muon track, excluding the muon itself, must be 
less than 10% of the muon p^. 

In the signal regions, the measurement of E!f ms is performed using all clusters of 
energy deposits in the calorimeter up to \rj\ of 4.5. The calibration of these clusters takes 
into account the different response of the calorimeters to hadrons compared to electrons 
or photons, as well as dead material and out-of-cluster energy losses [58, 59]. In the 
control regions, two additional definitions of E!f ms are used to account for the different 
treatment in the signal and control regions of electrons and muons. For control regions that 
involve electrons (W — > ei^+jets, Z — > e + e~+jets, explained in more detail below), E™ 1SS '^ 
is obtained by adding the electron clusters to the missing transverse momentum vector 
thereby removing the electron contribution to the calculation of E!f ms in order to match its 
definition in the signal regions where electrons are vetoed: E™ 1SS '^ = |p™ lss +p^ cctrons |. This 
yields missing transverse momentum which, as in invisible Z decays, does not take into 
account the decay products of the Z boson. The third E!f lss variant, E™ ss '^, is used in the 
exclusive W —> fiv+jets control regions and takes into account the muon contribution to 
E!f lss . It is the negative sum of the calorimeter-based E!f lss and the transverse momentum 
of muons, which do not deposit much energy in the calorimeters: E™ 1SS '^ = | p™ lss — p™ uons | . 
Without accounting for the muon contribution, the calorimeter-based E™ ss in W and Z 
decays involving muons mimics E™ ss in Z — > vv events. 

4 Event selection 

All data passing detector quality requirements are considered for the analysis. Events must 
be accepted by an inclusive E!f lss trigger [60, 61] that is found to be 98% efficient for events 
with EJ^ iss above 120 GeV, and more than 99% for Ef iss above 150 GeV. At 120 GeV, a 
small residual dependence on pile- up of the E™ 1SS trigger efficiency is found. Over the 
full 2011 dataset, where the pile-up varied from an average of 3 interactions per bunch 
crossing at the beginning of the year to 17 at the end of the year, an efficiency variation of 
1.5% is observed and a correction is applied to account for this variation. For E!f lss above 
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220 GeV, there is no measurable efficiency variation. Events are further required to satisfy 
a set of pre-selection and kinematic criteria that are aimed at selecting monojet events 
from good-quality pp collisions, as well as reducing electroweak, multijet, non-collision, 
and detector-induced backgrounds. These criteria require the event to have a monojet 
topology characterised by one unbalanced high-px jet resulting in large E™ 1SS : 

• A reconstructed primary vertex with at least two associated tracks (with > 
0.4 GeV) is required [62]. This ensures that the recorded event is consistent with 
a proton-proton collision. 

• The highest-pT jet must have a charge fraction / c h = YlPr* ^/Pr" > 0.02, where 
^^trackjet . g ^ e sca j ar sum Q f the transverse momenta of tracks associated with 
the primary vertex within a cone of radius R = 0.4 around the jet axis, and is 
the transverse momentum of the jet as determined from calorimeter measurements. 
Furthermore, events are rejected if they contain any jet with an electromagnetic 
fraction / cm (fraction of the jet energy measured in the electromagnetic calorimeter) 
of less than 0.1, or any jet in the pseudorapidity range \rj\ < 2 with / em > 0.95 and 
a charge fraction f c ^ < 0.05. These requirements suppress jets produced by cosmic 
rays or beam-background muons that interact in the hadronic calorimeter. 

• Additional selection criteria to reject events with significant detector noise and non- 
collision backgrounds are applied: events are rejected if any jet with p^ > 20 GeV 
and 1 7/| < 4.5 does not pass all of the additional quality criteria described in ref. [63]. 

• The leading jet has to be within \rj\ < 2, and no more than two jets with > 30 GeV 
and \rj\ < 4.5 are allowed. Back-to-back dijet events are suppressed by requiring the 
sub-leading jet not to point in the direction of p™ lss : | A0(p miss , p J ^ t2 )| > 0.5. 

• An electronics failure affecting 20% of the data sample created a small dead region in 
the second and third layers of the electromagnetic calorimeter. Any event with the 
two leading jets inside the affected region and either of the two jets pointing in the 
direction of E™ 1SS is removed from the sample to avoid fake signals. This condition 
removes only a few percent of the affected subset of the data. 

• Events are required to have no SR electron or muon. In the background control 
regions, electrons and muons are explicitly selected. The electron and muon selection 
criteria in the signal and control regions are given in section 3. 

Although the results of this analysis are interpreted in terms of the ADD model and 
WIMP pair production, the event selection criteria have not been tuned to maximise the 
sensitivity to any particular BSM scenario. To maintain sensitivity to a wide range of BSM 
models, four sets of overlapping kinematic selection criteria, designated as SRI to SR4, are 
defined (table 2). The requirement on the leading jet pt is the same as that on E!f lss 
for all signal regions. In comparison with the previous ATLAS monojet search [16], the 
veto on additional jets is less stringent, allowing a second jet in the event thereby reducing 
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Signal regions 


SRI SR2 SR3 SR4 


Common requirements 


Data quality + trigger + vertex + jet quality + 
\rji ctl \ < 2.0 + |A^»(p™ ss ,p^ t2 )| > 0.5 + iV jcts < 2 + 
lepton veto 


Ef iss , pJf 1 > 


120 GeV 220 GeV 350 GeV 500 GeV 



Table 2. Definition of the four overlapping signal regions SR1-SR4. 



systematic uncertainties from ISR/FSR (see section 6) and increasing signal selection effi- 
ciencies. The signal region with the lowest E™ 1SS requirement (SRI) is chosen such that the 
E™ 1SS trigger is nearly 100% efficient. The signal region with the highest E™ 1SS requirement 
(SR4) is chosen so that there remain enough events in data control samples to validate MC 
predictions and estimate SR backgrounds in a data-driven way. 

5 Background estimation 

A number of SM processes can pass the monojet kinematic selection criteria described 
above. These backgrounds include, in decreasing order of importance: Z and W boson plus 
jets production, single or pair production of top quarks, multijet production, cosmic-ray 
and beam-background muons 2 (collectively referred to as non- collision background [64]), 
and di-boson production (WW, WZ, ZZ). The dominant Z/W plus jets backgrounds are 
estimated using control regions in the data with corrections that account for differences 
between the selection requirements of the signal and control regions. The non-collision 
and multijet backgrounds are also estimated from data while the di-boson and top-quark 
backgrounds are obtained from MC simulations. 

5.1 Backgrounds from Z/W+jets 

The dominant background process for this search is irreducible and consists of the pro- 
duction of Z bosons in association with jets, where the Z decays to two neutrinos. A 
substantial source of reducible background is SM W boson plus jets production where 
the W decays to a charged lepton (r, e, or ^ in decreasing order of importance) and a 
neutrino. This process leads to a monojet final state if the lepton is outside the detector 
acceptance or is missed because of reconstruction inefficiencies. The Z and W boson plus 
jets backgrounds, collectively referred to in the following as electroweak backgrounds, are 
determined in a data-driven way: 

1. Control regions are defined by explicitly selecting electrons or muons while keeping 
the same jet and E™ 1SS selection criteria as in the signal regions. Corrections are 
made for contamination of these control regions from processes other than Z or W 
decays. Note that the irreducible background from Z — > z^z/+jets must be treated 

2 Originating either from protons going in the direction of the experiment and hitting the LHC collimation 
system or gas molecules in the beam-pipe near the ATLAS interaction point. 
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SR process 


W — > r^+jets Z — > t + t +jets 
Z -> vu+jets W -> ev+iets + _ . 

W — y /iz^+jets Z — y fi , u +jets 


CR process 


W — > eu+jets 

W — > [iv+jets ,„ , „ 

, W — > /ii/+jets W — > e^+jets Z — > fj, r fj, +jets 
Z — > e + e +jets 

Z — > /i + ^~+jets 



Table 3. Overview of processes in the control regions (CR) used to estimate background contribu- 
tions to processes in the signal regions (SR). 



specially: E™ 1SS is redefined in the CRs that involve electrons and muons, in order to 
mimic the kinematics of the signal regions (see below). 

2. Corrections are applied to account for differences in trigger and kinematic selection 
criteria between the control and signal regions. The control-to-signal region transfer 
factors are obtained using both data and simulation. 

In this approach, the modelling of the jet and E!f lss kinematics of the electroweak back- 
grounds is obtained directly from data. Simulations are therefore used only for quantities 
related to the electron and muon selection criteria, and only through ratios where sys- 
tematic uncertainties related to the jet and E™ 1SS selection criteria of the control regions 
cancel. Theoretical uncertainties normally associated with MC estimates are significantly 
reduced; only distributions related to the electron and muon selection criteria have to be 
well modelled in the simulations. Further experimental uncertainties that impact the back- 
ground prediction, such as the jet energy scale (JES) and resolution (JER) [55], the trigger 
efficiency, and the luminosity measurement [36, 37], are minimised by this approach. 

The control regions are expected to have no contamination from BSM signals that 
would normally pass the monojet event selection criteria. They are chosen such that 
they are dominated by Z and W decays with reconstructed electrons or muons. The 
selection criteria follow closely those used in Z and W cross-section measurements [65] . The 
kinematic selection criteria on E!f lss and jet px of the signal regions are also applied. As a 
result, each CR is split into a subset of four regions corresponding to the four signal regions. 
Four visible decay modes contribute to the control regions: W — > ez^+jets, W — > /uz^+jets, 
Z —7- e + e _ +jets, Z — > ^ + /i~+jets. Using these, all contributions to the signal regions from 
Z and W decay modes are estimated with the same method (except for Z — > e + e _ +jets, 
which is found to be negligible in the signal regions because both e + and e~ would have to 
be missed in the event selection). In total, six background processes in each signal region 
are predicted based on four control region processes, as detailed in table 3. For control 
regions that include processes with electrons, an electron trigger is used that requires a 
correction to account for differences in efficiency and acceptance compared to the E™ 1SS 
trigger used for the signal regions. Note that the acceptance is defined as the ratio of the 
number of events within the detector volume that pass analysis requirements to the number 
of originally simulated events. The efficiency is defined as the ratio of the number of events 
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within the detector volume at reconstruction level to that at the original simulation level. 
A different trigger is required for control regions that include electron selection criteria 
because, for the inclusive E!f lss trigger, the energy deposited by electrons is included in 
the E™ 1SS measurement at trigger level and results in the selection of a different kinematic 
region than that of the signal regions, which exclude electrons. Muons, however, do not 
deposit large amounts of energy in the calorimeters and are not explicitly included in the 
£;miss trigger. The specific selection criteria for the four control region processes are given 
in the following: 

• W — > ez^+jets: Events are selected using electron triggers with thresholds of 20 or 
22 GeV depending on the data-taking period. The CR electron definition is used (see 
section 3) and exactly one electron with a px of at least 25 GeV is required. Events 
with additional electrons or muons are discarded. All triggers used are fully efficient 
above the chosen px cut value. If an object is reconstructed as both an electron 
and a jet, the jet is removed from the reconstructed jet collection if Ai?(e, jet) < 0.2 
while the electron is kept. To further improve the W purity, E™ 1SS > 25 GeV and 
40 < m,T < 100 GeV are required, mx is the transverse mass and it is defined 
as rriT = \J 2pp E™ ss (1 — cos A(j)(£, E™ 1SS )), using the px of the lepton (electron or 
muon). As mentioned earlier, the same selection criteria on jet px and E!f lss are 
applied in the control regions as in the signal regions (see table 2). However, when 
the W — > ez/+jets CR is used to estimate the contribution of Z — > vv+jets to each SR, 
a special CR is defined where E™ ms is substituted by S™ 1SS '^ to mimic the kinematics 
of the decay of the Z boson to two undetected neutrinos. The standard calorimeter- 
based E™ ms is used for the CR to estimate the W — > e^+jets contribution to the 
SRs. 

• W — > /x^+jets: Events have to pass the same inclusive E™ 1SS trigger that is used for 
the signal regions. Exactly one CR muon (see section 3) is required, and events with 
additional electrons or muons are rejected. Cuts on the transverse mass and missing 
transverse momentum are applied to improve the purity for W's: mx > 40 GeV, 
^miss,/i ^ 25 Q e y Note that the E!f lss that includes the muon contribution, E™ ss '^, 
is used for the PF-specific selection cuts. For each kinematic region listed in table 2, 
the standard calorimeter-based E!f lss is used to define the CRs for the estimates of 
both Z — > z/P+jets and W — > pi/+jets in the corresponding SR. 

• Z — > e + e _ +jets: Electron triggers are used in this channel. Exactly two opposite-sign 
electrons are required and events with additional electrons or muons are discarded. 
The selected electrons have to satisfy pr > 25 (20) GeV for the leading (sub-leading) 
electron. Jet-electron overlap removal is performed as described above for W —> 
ez^+jets. Finally, to enhance the fraction of Z's, an invariant mass requirement of 
66 < m e + e - < 116 GeV is applied. The modified missing transverse momentum 
gmmsrf . g uge( j ^ ,-ief^e these CRs, which are used to estimate the Z — > vv+jets 
contribution to the SRs. 
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• Z — > yU + /i _ +jets: The inclusive i?™ lss trigger is used in this channel. Exactly two 
opposite-sign CR muons (defined in section 3) are required and events with additional 
electrons or muons are rejected. An invariant mass of 66 < rn^+^- < 116 GeV is re- 
quired to select Z candidates. The signal-region selection criteria are then applied on 
the calorimeter-based £?j? lss , for the Z — > i/P+jets, Z — > T + T _ +jets, Z — > /i + /i~+jets 
estimates. 

Using the control regions defined above, the background contribution to the signal regions 
for each combination of CR and SR processes mentioned in Table 3 is derived using: 



^predicted = (Ar Data _ ^Bkgj . (j . j, 



)6 PI f; ss ' rpmiss 7\rMC 
X ~ A trig r X l^MC • ^ A > 

Ac e r e z/w -e t ■ L t iV je t/E" 



Data in the control regions are corrected for contamination arising from other sources. 
Correction factors (C) based on MC simulation and data are applied together with the 
transfer factor T to obtain the number of background events, A^ cdlctcd , predicted in the 
signal region. The terms appearing in the second line of equation 5.1 are: 

• iV(j^ ta and A r ^ ,ltl ^ et are the number of data and multijet events in the control region, 
respectively. To estimate the multijet contamination of control regions by processes 
with identified electrons, the selection cuts, in particular the isolation cuts, are var- 
ied. The fake rate in those regions is extracted from data using real and fake electron 
efficiencies determined from samples enriched in electrons and jets. Using this esti- 
mate, the multijet contamination is predicted to account for 1-2% of the events in the 
W — > ez/+jets electron control region when predicting the Z — > vv+jets contribution 
to the SRs. For other control regions containing electrons or muons, the multijet 
contamination is found to be negligible using similar techniques. 

• /ew is the estimated fraction of events, after multijet corrections, due to contami- 
nation of the control region by other electroweak or other SM processes. This con- 
tamination is due to top-quark and di-boson decays as well as decays of Z or W 
bosons to leptons of a flavour other than the one selected for that control region. 
The contribution of this contamination is obtained from MC simulation and is about 
2% for Z bosons, and 10% for W bosons. The top-quark and di-boson contribution 
is negligible. As explained above, using ratios of MC estimates (/ew i n this case) is 
advantageous as it leads to cancellations of systematic uncertainties. 

• Ai and are the lepton acceptance obtained from simulation and the identification 
efficiency obtained from data, respectively. 

• ez/w are the efficiencies for the Z or W boson selection criteria obtained from sim- 
ulation. The factors A%, eg, and ez/w correct for the fact that leptons and Z /W 
bosons are required only in the control regions. 
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• e^ ng and Li are the electron trigger efficiency (obtained from data) and the corre- 
sponding luminosity associated with this trigger for the relevant control region. For 
muon control regions these factors do not apply because the signal-region trigger is 
used in the definition of the CR. 

• e ^ i m and L E mi SS are the E!f lss trigger efficiency (obtained from data) and the cor- 
responding luminosity, and are only relevant for electron control regions where the 
electron trigger efficiency and luminosity (e/ lg and Lg) need to be corrected, account- 
ing for the different triggers used in the definition of the signal and control regions. 

7V MC 

• The transfer factor T = j^mc 11 — is the ratio of simulated background events in 

the signal region (for example Z — > vv+jets) to simulated events of a control-region 
process (for example W — > ez^+jets) with only jet and E™ ss related selection require- 
ments applied. This term translates the number of observed events in the CR in the 
data to the predicted number of events in the signal region. Depending on the control 
region and the signal-region background component being determined, this factor can 
account for ratios of branching fractions, ratios of W+jets to Z+jets cross sections, 
and phase-space differences between the control and signal regions for a given source 
of background. 

The correction factors and electroweak background predictions in equation 5.1 are deter- 
mined in bins of E™ 1SS for the final background prediction, and in bins of the leading and 
sub-leading jet pt for the jet-px plots in figure 2. 

5.2 Multijet backgrounds 

Multijet events where one or more jets are severely mismeasured constitute a background 
that is not well modelled in the simulation. In order to measure this background from 
data, a sample is selected by applying all signal-region selection criteria except for the jet 
vetoes: either a second jet with | A(j)(Erf lss , p J ^ t2 )| < 0.5 is required, or the third-jet veto is 
reversed by requiring A?~ jet = 3, | A(j)(Ef iss , p J ^ t3 )| < 0.5 and | A(j)(Ef iss , p J ^ t2 )| > 0.5. These 
two samples are used to predict the multijet background from the resulting di- or trijet 
events. Contributions to these event samples from top-quark, and Z or W production are 
subtracted. The MC simulation is used for the top-quark contribution. For Z and W, 
MC estimates normalised to data are used for the subtraction. The multijet background is 
then estimated by fitting a straight line to the second or third jet distributions in events 
passing the two selection criteria and then extrapolating the fit below a pi of 30 GeV. For 
this value of the transverse momentum, the jets fall below the threshold and can pass the 
monojet selection criteria. Note that the number of trijet events where both sub-leading 
jets are mismeasured, and fall below the jet threshold, is negligible compared to the case 
where either the second or the third jet is lost. The resulting background estimates are 
given in table 5. They are at most 1% of the total background predicted for SR1-SR3, and 
are negligible for SR4. 
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5.3 Non-collision backgrounds 

Non-collision backgrounds in the signal regions are estimated using a dedicated algorithm 
that identifies beam-background muons that go through the detector along the direction 
of the beams. The algorithm selects through-going muons based on timing information 
obtained from the muon chambers in the forward regions. It combines this information 
with calorimeter energy clusters by matching them in (j). Unpaired proton bunches, where 
the bunch from one of the proton beams is empty, are used to determine identification 
efficiencies of the algorithm for beam-background and cosmic-ray muons. These efficien- 
cies (e^ag~ CoU )> typically 20-50%, are used together with the number of beam-background 
{halo) candidates found in the signal regions, to predict the level of non-collision back- 
ground (iV non _ co n = A^aioAtag _CoU )■ More details of this background component are 
given in ref. [64]. It contributes mainly in SRI and SR2 at less than 1%. The predictions 
are given in table 5. 

5.4 Systematic uncertainties on background estimates 

The dominant systematic uncertainties associated with the electroweak background es- 
timates are on JES and E™ ss , as well as theoretical uncertainties on the shape of W 
kinematic distributions and the ratio of Z and W plus jets production cross sections. The 
latter theoretical uncertainty is relevant because background predictions from W control 
regions are also used to estimate Z —> vv+]ets contributions to the signal region. Addi- 
tional systematic uncertainties are due to the muon momentum scale and resolution, the 
data-driven scale factors for lepton trigger and reconstruction efficiencies, statistical uncer- 
tainties associated with the limited size of MC samples, the subtraction of the electroweak 
contamination (/ew)i and, in the electron control regions, the multijet contamination. 
The uncertainties from pile-up variations are found to be negligible, as are those from the 
electron energy scale, resolution, and JER. 

The systematic uncertainties associated with the small multijet background are esti- 
mated to be 100%. They are obtained by changing the fit range for the pt extrapolation 
and varying the scale factors for the Z and W background prediction by 10%. All variations 
are within a factor of two of the central predictions. 

Systematic uncertainties on the non-collision background are 10%. This estimate cor- 
responds to the average fraction of unpaired proton bunches that are used to determine 
e"°g _co11 ', and that are close (separated by 25 ns in time) to an unpaired bunch in the 
opposite beam. Such configurations may lead to double counting in the efficiency estimate, 
and their total contribution is hence considered as an uncertainty. 

The JES and JER uncertainties are evaluated using a combination of data-driven and 
MC-based techniques [55]. These methods account for the variation of the uncertainty 
with jet pt and r], and for the presence of nearby jets. The E!f lss uncertainty is derived 
from the JES and JER uncertainties by propagating the relative jet-level variations to 
the calorimeter cluster based E™ 1SS . Since ratios are used to extrapolate from the control 
regions to the signal regions, the effects of these uncertainties tend to cancel. 
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Table 4. Relative systematic uncertainties for all signal regions (in percent). Individual contri- 
butions are summed in quadrature to derive the total numbers. The MC statistical uncertainty is 
included in the total systematic uncertainty. 

Theoretical uncertainties on Z and W production and the shape of W kinematic dis- 
tributions are evaluated by comparing background estimates using kinematic Z/W distri- 
butions from different generators (ALPGEN and SHERPA). Uncertainties on /ew are derived 
by comparing ALPGEN to PYTHIA [66], but also by taking into account JES and lepton-scale 
uncertainties. 

Systematic and statistical uncertainties on all background estimates are given in ta- 
ble 4. The contribution from lepton scale factors is the quadratic sum of electron and muon 
uncertainties. The uncertainties from di-boson, top-quark, multijet, and non-collision back- 
grounds are summed in quadrature. A 20% uncertainty is assigned for the di-boson and 
top-quark MC-based estimates. This value is dominated by the JES uncertainty (16%), 
but also takes into account uncertainties of the trigger efficiency, luminosity measurement, 
and lepton identification uncertainties. 

5.5 Background summary and additional checks 

An overview of all backgrounds is given in table 5 (cf. table 3 for the definition of the 
control regions). The final Z — > uv+}ets predictions are estimated from a combination of 
the predictions of the four control regions. The combination is the error- weighted average 
calculated taking into account correlations of uncertainties. The combined Z — > z^+jets 
prediction is dominated by the W control-region estimates and based on the assumption 
that the ratio of Z+jets to VK+jets cross sections is well modelled in the simulation. This 
assumption is supported by dedicated measurements [66] , albeit for smaller jet momenta 
than the ones used in SR2 to SR4. The theoretical uncertainty on the ratio oiZtoW cross 
sections is included in the uncertainty derived from comparisons of different MC generators 
discussed above. 

The electroweak background estimate, which relies on an exclusive W or Z selection 
in the control regions, is compared to two alternative correlated methods. In the first of 
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Table 5. Overview of predicted SM background and observed events in data for 4.7 lb -1 for 
each of the four signal regions. The total uncertainty quoted is the quadratic sum of statistical 
and systematic uncertainties. Observed and expected 90% and 95% CL upper limits on the non- 
SM contribution to all signal regions are also given in terms of limits on visible cross sections 
{a vis = a x A x e). The 90% CL upper limits are given to facilitate comparisons with other 
experiments. 

these, which was the main method used in the previous ATLAS monojet search [16], an 
inclusive control region is defined by only inverting the lepton veto while keeping all other 
selection criteria the same as in the signal regions. No additional Z- or W^-specific invariant 
or transverse mass selection criteria are applied, thereby yielding a mixed control sample 
dominated by W and Z bosons. The resulting background predictions are found to be 
consistent with those of the default method. The second alternative modifies the lepton 
definition in the control regions. Instead of applying lepton selection criteria in control 
regions that are more stringent than those of the signal regions, a modified exclusive control 
region is defined. The selection criteria include less stringent lepton definitions where the 
lepton veto cuts of the signal region are simply inverted, and dedicated Z or W selection 
criteria are used. These background predictions are also found to be consistent with the 
default method. 

Distributions from all four visible decay modes used to determine the background in 
SRI (labeled CR1) are shown in figure 1. The distributions are obtained by applying the 
exclusive Z and W selection criteria plus SRI kinematic cuts on E!f lss and jets, as well 
as vetoes on additional electrons or muons. It should be noted that shape differences in 
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the E™ 1SS distributions between data and MC are irrelevant for an accurate background 
prediction in the signal regions, because the E!f lss distribution obtained from control-region 
data is used directly to predict the backgrounds in the signal regions. Distributions of 
variables that are subject to MC-based efficiency or acceptance corrections, namely those 
involving electrons or muons, need to agree in shape between data and MC (see figure 1, 
where good shape agreement is found for the leading electron and muon distributions). 
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Figure 1. Kinematic distributions in the control regions corresponding to SRI (labelled CR1) are 
shown. The upper row is the leading electron and muon px distribution for Z — > e + e~+jcts (left) 
and Z — > /i + /i~+jcts (right) and shows distributions after SRI cuts on jets and E^ lss . The lower 
row is the missing transverse momentum distribution E'™ ss ' 9 ' for W — > e^+jets (left) and i?™ lss for 
W -t /w+jcts (right) also after SRI jet and E™ lss cuts. 
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6 Results and interpretation 
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Figure 2. Kinematic distributions for signal regions SRI on the left and SR4 on the right. Signal 
distributions for ADD and WIMP samples for cross sections equal to the excluded values are drawn 
as dashed lines on top of the predicted background distributions. The electroweak backgrounds 
(see equation 5.1) are determined in bins of the variable that is plotted. 



The SM predictions are found to be consistent with the number of observed events 
in data for all signal regions considered. Comparisons of the SM predictions to the mea- 
sured E™ 1SS and leading and sub-leading jet pr distributions are shown for SRI and SR4 
in figure 2. For illustration, the figures also contain simulated signal distributions for ADD 



- 18 - 





Dlxl [ /o J 


odo r oy i 
orxz [ /o J 


orto [ /o J 


Ol\,4 [ /0 J 


Z — y z/z^±jets 


1. / UD It U.UlO 


n i en _l n nn/i 
u.ioy it u.uu4 


n m 7n 4- n nm 

U.U1/U It U.UUlO 


O 0097 4- O OOO^ 
U.UUZl It U.UUUO 


A F>F> n 9 

A1J1J, 71 — Z 


Qfi q -i- n 9 
ou.y it u.z 


Q 9 4- n 1 

y.z ± u.i 


9 fin 4- n 07 
Z.DU It u.u / 


O 7/1 4- O O/l 
U. 1 4 It U.U4 


A T^T^ m Q 

AJJJJ, 71 — o 


OO.Z ± U.Z 


11 7 4_ n 1 

11. 1 ± U.J. 


■? 09 4- n oa 
o.yz ± u.uo 


1 18 4-0 O^ 
J.Jo it U.UO 


A PiT^ m /I 

AJJJJ, 71 — 4 


04.0 ± U.Z 


lo.o it U. 1 


A 07 4- O OQ 

4.y / ± u.uy 


1 fi7 4- O O^ 
1.0/ It U.UO 


IWJL), 71 — 


oO.l ± u.z 


i/i ^ 4- n 1 

14.0 It U. 1 


c;n 4- O OQ 

o.ou ± u.uy 


9 OO 4- O Ofi 
Z.UU It U.UD 


A in fi 

IWJL), 71 — 


oO.U ± U.Z 


10. u ± u.z 


fi Ol 4-010 
D.Ul it U.1U 


9 90. 4- O Ofi 
Z.Zo It U.UD 


D 1 . 771-^. — 1U Kjev 


9n c-i-in 
zu.O ± u.o 


o qin i 
o.o it U.I 


O KA 4- O Ol 
U.04 it U.U1 


O OQ 4- O 01 

u.uy it u.ui 


L/l, 774-^ — 1UUU <orev 


OZ.Z It U.4 


i n t. 4- n 9 

1U.O It u.z 


9 88 4- O O/l 
Z.oo It U.U4 


O 7Q 4- O 09 

u. ( y it u.uz 


pic: 1 P oV 

JJD, 771-^ — 1U bcv 


Qfi /i 4- n a 

OU.4 It U.4 


O.O It U.Z 


9 O/l 4- O OT. 
Z.U4 It U.UO 


U.OZ It U.UI 


n5 m — 1000 OeV 




12fi + n2 

X & • I V./.— j 


4 1 4 + 05 


i 94 4- o fn 


D9, m x = 10 GeV 


36.9 it 0.5 


12.9 ib 0.3 


4.23 ±0.15 


1.31 ±0.08 


D9, 77i x = 1000 GeV 


37.6 it 0.5 


13.9 ±0.3 


4.70 ±0.16 


1.68 ± 0.09 


Dll, m x = 10 GeV 


30.3 ± 0.4 


12.3 ±0.3 


4.57 ±0.15 


1.52 ±0.09 


Dll, m x = 1000 GeV 


33.7 ±0.5 


17.0 ±0.3 


7.56 ± 0.20 


3.27 ±0.13 



Table 6. Typical acceptances determined with MC simulations for the main background process 
Z — > i/jz+jets as well as for ADD and selected WIMP samples. For the Z — > w+jets sample at 
least one parton with a minimum transverse momentum of 20 GeV is required, for the ADD and 
WIMP samples it is at least one parton with a momentum of 80 GeV. The values are given in 
percent and errors are statistical only. 



and WIMP models added to the total background. Agreement both in the shape and 
the overall normalisation between SM predictions and data is observed in all cases. To 
facilitate comparisons with other experiments both 90% and the more conventional 95% 
confidence level (CL) upper limits are produced. These limits are on the visible cross sec- 
tion defined as cross section times acceptance and efficiency (a x A x e) and they are based 
on the modified frequentist CL S prescription [67] . The limits are derived by comparing the 
probabilities, based on Poisson distributions, that the observed number of events is com- 
patible with the SM and the SM-plus-signal expectations. The mean values of the Poisson 
distributions are determined by the signal prediction, plus contributions from background 
processes extrapolated from the CRs to each SR. The number of events is integrated over 
the whole SR. Expected limits are obtained by repeating the analysis with pseudo-data 
obtained from Monte Carlo simulations. The distributions of the simulated probabilities 
for many pseudo-experiments allow ±1<t bands to be plotted for the expected values. Sys- 
tematic uncertainties (and their correlations) associated with SM backgrounds and the 
integrated luminosity are taken into account via nuisance parameters using a profile like- 
lihood technique [68]. The nuisance parameters are assumed to be Gaussian distributed 
and are constrained in the fit to lie within ±1<t of their value determined in the study 
of the systematic uncertainties. The resulting visible cross-section limits, which apply for 
any source of BSM events, are summarised in table 5. Typical efficiencies of selection cri- 
teria related to jets and E^ ms of e ~ 83% are found in simulated Z — > z^P±jets, WIMP 
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or ADD samples. Typical acceptances are given in table 6. The negative search results 
are interpreted in terms of limits on ADD and WIMP model parameters in the following 
sections. 

6.1 Large extra dimensions 




M D [GeV] Number of extra dimensions 

Figure 3. Left: Visible cross sections in SR4 as a function of Md as predicted by the effective ADD 
theory, for n = 2, 4, 6 extra dimensions. The coloured bands correspond to the theoretical systematic 
uncertainties (PDF, ISR/FSR, scale). The horizontal lines are the expected and observed cross- 
section limits at 95% CL, taking into account experimental systematic uncertainties fully correlated 
between signal and background, as well as uncertainties on the luminosity estimate, trigger efficiency, 
and MC statistical uncertainties. Right: 95% CL lower limits on Md for different numbers of extra 
dimensions based on SR4. Observed and expected limits including all but the theoretical signal 
uncertainties are shown as solid and dashed lines, respectively. The grey ±f<r band around the 
expected limit is the variation expected from statistical fluctuations and experimental systematic 
uncertainties on SM and signal processes. The impact of the theoretical uncertainties is shown by 
the red small-dashed ±lcx limits. The previous ATLAS limit [16] is also shown for comparison. 

Experimental and theoretical systematic uncertainties that affect the ADD signal are 
considered in order to set limits on the model parameters. The experimental uncertainties 
on JES, JER, and E™ 1SS are considered to be fully correlated with those obtained for 
the background estimate. They range from 3-10% depending on the signal region and 
the number of extra dimensions. An additional 1% uncertainty on the trigger, and a 
3.9% uncertainty on the luminosity are also considered for the signal simulation only. 
Theoretical uncertainties on the expected ADD signal are associated with the PDF set, 
ISR/FSR, and the factorisation and renormalisation scales. For the PDF uncertainties, 
the CTEQ6 . 6 error sets are used, converted from 90% to 68% CL. They range from 4- 
14% on the product of signal cross section and acceptance (a x A), depending on the 
number of extra dimensions. Uncertainties coming from the modelling of ISR/FSR are 
determined by varying the simulation parameters of PYTHIA within a range that is consistent 
with experimental data [69]. The resulting uncertainties vary from about 3-14%. The 
dominant theoretical systematic uncertainty affecting mostly the cross section rather than 
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Table 7. 95% CL lower (upper) limits on Md (R) for n=2-6 extra dimensions, using a dataset 
corresponding to 4.7 fb _1 at = 7 TeV. These results are obtained using the selection criteria 
of SR4. All values correspond to the nominal observed limits excluding theoretical uncertainties 
in figure 3. The last two columns show the relative difference between the full cross sections and 
those of the truncated phase space (s < M^). The ADD cross sections are calculated at both LO 
and NLO, and the limits are derived from the full, not the truncated, phase space. 



the acceptance is from the factorisation and renormalisation scales. Varying these scales 
between twice and half their default values results in 20-30% uncertainties on a x A. 3 

The visible cross sections predicted by the ADD generator for SR4 are shown for 
n = 2, 4, 6 extra dimensions as a function of Md on the left-hand side of figure 3. Theoretical 
systematic uncertainties are shown as coloured bands around the cross-section curves. The 
95% CL expected and observed limits on the visible cross section a x A x e are shown as 
horizontal lines. The effect of restricting the simulated phase space to the kinematic region 
where the ADD effective field-theory implementation is valid is probed by evaluating the 
cross section after discarding events for which the parton centre-of-mass energy s > M£. 
The amount by which the truncated cross sections differ from the full ones provides a 
measure for the reliability of the effective field theory. This difference increases from SRI 
to SR4 and with the number of extra dimensions. While the model with n = 2 extra 
dimensions is found to be insensitive to truncation effects for Md values near the resulting 
limits for all signal regions, n = 3, 4, 5, and 6 extra dimensions show differences of 2%, 10%, 
30%, and 50% between full and truncated cross sections for SR4 and Md values close to the 
actual limits (table 7). This demonstrates that the high energy and integrated luminosity 
used in this search allow to probe kinematic regions where the effective field-theory model 
is not entirely valid. 

The 95% CL lower limits on Md versus n for the full phase space, not the truncated 
one, are shown for SR4 on the right-hand side of figure 3. The selection criteria of SR4 
provide the best expected limits and are therefore used here. Limits from SRI, SR2, SR3 
are typically 35%, 15%, 5% worse, respectively. The expected and observed limits in figure 3 
are produced taking all but the theoretical uncertainties into account. The grey ±lc band 
around the expected limit shows the variation anticipated from statistical fluctuations 
and from experimental systematic uncertainties on background and signal processes. The 

3 Note that in ref. [16] the squared factorisation and renormalisation scales were varied between twice 
and half their default values. 
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impact of the theoretical uncertainties associated with PDFs, ISR/FSR, and factorisation 
and renormalisation scales is represented in the right-hand panel by dashed ±lcr lines on 
either side of the observed limit. The resulting limit is taken as the observed line excluding 
theoretical uncertainties. 4 All limits from SR4 are summarised in table 7, where the lower 
(upper) limits on Md (R) are shown for cross sections calculated at LO and NLO. The 
K-factors (defined as onlo/clo) f° r w = 2, 3, 4, 5, 6 extra dimensions are 1.20, 1.20, 
1.17, 1.13, 1.09, respectively, and have been derived for the selection criteria of SR4 by the 
authors of ref. [49]. M D values below 4.17 (4.37) TeV for n = 2 and 2.51 (2.53) TeV for 
n = 6 are excluded at 95% CL at LO (NLO). 

6.2 WIMP pair production 

Systematic uncertainties on WIMP pair production are treated similarly to those of the 
ADD limits, except for the PDF and ISR/FSR uncertainties. The former are determined 
using CTEQ6M error sets for the relative uncertainty around the CTEQ6L1 central value. The 
ISR/FSR uncertainties are estimated differently in a way that is appropriate for the high- 
Pt ISR/FSR regime probed here: a WIMP pair recoils against a high-pr ISR/FSR jet, 
whereas for ADD, additional low-px ISR/FSR jets dominate the uncertainty due to the 
impact of the jet veto. 

The JES/JER/£^™ 1SS experimental uncertainties lead to 1-20% uncertainties on the 
WIMP event yield depending on the signal region and the effective operator considered. 
Other experimental uncertainties affecting the WIMP event yield are associated with the 
trigger efficiency (1%) and the luminosity measurement (3.9%). The ISR/FSR uncertain- 
ties are estimated by varying the jet matching scale between MADGRAPH5 and PYTHIA by 
a factor of one half and two. Moreover, the a s scale is varied in PYTHIA within a range 
that is consistent with experimental data [69]. The resulting uncertainties on a x A, added 
in quadrature, range from 3-5% for the matching scale and 4-6% for a s depending on 
the signal region. A negligible dependence of the ISR/FSR uncertainties on the choice of 
effective operator is found. PDF uncertainties impact mostly the signal cross section and 
hardly the acceptance. They are found to depend on the effective operator chosen and not 
the particular signal region (since overall cross-section differences affect the signal regions 
in the same way). Uncertainties ranging from 4% and 5% for operators D9 and D5 to 
16% and 18% for Dll and Dl are found. As for the ADD model, the dominating theoret- 
ical systematic uncertainty is from the factorisation and renormalisation scales. Varying 
these scales between twice and half their default value results in 30% signal uncertainties, 
independent of the effective operator choice or the signal region. 

Figure 4 shows the 90% CL lower limits on the suppression scale M*, for all operators 
probed as a function of WIMP mass m x . These limits on M* are derived from the cross- 
section limits at a given mass m x . The values displayed are for the signal regions with 
the best expected limits, where those limits from SR3 and SR4 are typically within a few 
percent of each other, and those from SR2 (SRI) are 15-20% (40-50%) smaller than in SR3 

4 The previous ATLAS monojet search [16] has determined ADD parameter limits in a slightly different 
way. The effect of the signal cross section theoretical uncertainty was folded into the quoted limit and was 
not shown separately. 
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or SR4. The lower limits are based on simulation samples produced for m x between 10 and 
1300 GeV. Extrapolations are shown down to m x = 1 GeV. These are valid (and could be 
continued as constants to even smaller m x values entering the warm or hot dark-matter 
regime) since there is negligible change in cross section or kinematic distributions at the 
LHC for low-mass WIMPs. As before, the central values of observed and expected limits 
on M* are displayed taking into account experimental but not theoretical uncertainties. 
The effect of ztla variations on the expected limit due to statistical fluctuations and exper- 
imental uncertainties is shown as a grey band. The impact of the theoretical uncertainties 
is represented by dotted red ±l<r lines on either side of the observed limit. The nominal 
observed limit line excluding theoretical uncertainties is the final result. All values of the 
lower limits on the suppression scale M* at 90% and 95% CL are listed in table 8. For all 
operators, the lower limits are fiat up to m x = 100 GeV and worsen around m x = 200 GeV. 
Note that the M* limits for Dl are much smaller due to the inclusion of a factor m q /M* 
in the definition of the operator (see table 1). 
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Figure 4. ATLAS lower limits at 90% CL on M* for different masses of \ — the region below 
the limit lines is excluded. The 90% instead of the 95% CL lower limits are plotted because the 
former are used in the following figures 5 and 6. Observed and expected limits including all but 
the theoretical signal uncertainties are shown as dashed black and red solid lines, respectively. The 
grey ±1<t band around the expected limit is the variation expected from statistical fluctuations and 
experimental systematic uncertainties on SM and signal processes. The impact of the theoretical 
uncertainties is shown by the thin red dotted ±lrj limit lines around the observed limit. The M* 
values at which WIMPs of a given mass would result in the required relic abundance are shown as 
rising green lines (taken from [32]), assuming annihilation in the early universe proceeded exclusively 
via the given operator. The shaded light-grey regions in the bottom right corners indicate where 
the effective field theory approach breaks down [32]. The plots for Dl, D5, D8 are based on SR3, 
those for D9 and Dll on SR4. 
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Table 8. ATLAS 90% (95%) CL observed lower limits on the suppression scale M* as a function 
of WIMP mass m x . All values are given in GeV and correspond to the nominal observed limit 
excluding theoretical uncertainties. The signal regions with the best expected limits arc quoted in 
all cases, SR3 is used for Dl, D5 and D8, SR4 for D9 and Dll. 

The light-grey shaded regions in figure 4 indicate where the effective field theory ap- 
proach for WIMP pair production breaks down [32] (bottom-right corner in all plots). 5 
Except for some of the m x = 1300 GeV points, the M* limits set in this analysis are well 
above these bounds. No further measures are taken to ensure that the energy transfer in 
monojet events in this dataset remains in the valid region of the effective field theory. Such 
a region of validity cannot be defined without precise knowledge of the BSM physics, over 
which the effective operators integrate. 

Figure 4 also includes thermal relic lines (taken from [32]) which correspond to a 
coupling, set by M*, of WIMPs to quarks or gluons such that WIMPs have the correct 
relic abundance as measured by the WMAP satellite [30], in the absence of any other 
interaction than the one considered. Under the assumption that DM is entirely composed 
of thermal relics, ATLAS limits on M* that are above the value required for the thermal 
relic density exclude the case where DM annihilates exclusively to SM particles via the 
corresponding operator. Should thermal relic WIMPs exist in these regions (above the 
thermal relic line), there would have to be other annihilation channels or annihilation via 
other operators in order to be consistent with the WMAP measurements. 

In the effective operator approach, the ATLAS bounds on M* for a given m x can 
be converted to bounds on WIMP-nucleon scattering cross sections, which are probed by 
direct dark matter detection experiments. These bounds describe scattering of WIMPs 
from nucleons at a very low momentum transfer of the order of a keV. Depending on the 
type of interaction, contributions to spin- dependent or spin-independent WIMP-nucleon 
interactions are expected. The translation of ATLAS limits to bounds on WIMP-nucleon 

a Compared to ref. [32] the valid region of Dl shown here accounts for the factor of m q in the definition 
of Dl (see table 1). 
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Figure 5. Inferred 90% CL ATLAS limits on spin- independent WIMP-nucleon scattering. Cross 
sections are shown versus WIMP mass m x . In all cases the thick solid lines are the observed 
limits excluding theoretical uncertainties; the observed limits corresponding to the WIMP-parton 
cross section obtained from the — lotheory l mes m figure 4 are shown as thin dotted lines. The 
latter limits are conservative because they also include theoretical uncertainties. The ATLAS limits 
for operators involving quarks are for the four light flavours assuming equal coupling strengths 
for all quark flavours to the WIMPs. For comparison, 90% CL limits from the XENON100 [70], 
CDMSII [71], CoGcNT [72], CDF [19], and CMS [21] experiments are shown. 



scattering cross sections is done using equations (3) to (6) of ref. [32], and the results are 
shown in figures 5 and 6. 6 As in ref. [32] uncertainties on hadronic matrix elements are 
neglected here. The spin-independent ATLAS limits in figure 5 are particularly relevant in 
the low m x region (< 10 GeV) where the XENON100 [70], CDMSII [71] or CoGeNT [72] 
limits suffer from a kinematic suppression. Should DM particles couple exclusively to 
gluons via Dll, the collider limits would be competitive up to m x of about 20 GeV, and 
remain important over almost the full m x range covered. The spin-dependent limits in 
figure 6 are based on D8 and D9, where for D8 the M* limits are calculated using the D5 
acceptances (as they are identical) together with D8 production cross sections. Both the 
D8 and D9 cross-section limits are significantly smaller than those from direct-detection 
experiments. 

As in figure 4, the collider limits can be interpreted in terms of the relic abundance 

6 There is a typographical error in equation (5) of ref. [32] (cross sections for D8 and D9). Instead of 
9.18 x lCT 40 cm 2 the pre-factor should be 4.7 x 10" 39 cm 2 . 
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Figure 6. Inferred 90% CL ATLAS limits on spin-dependent WIMP-nucleon scattering. Cross 
sections are shown versus WIMP mass m x . In all cases the thick solid lines are the observed limits 
excluding theoretical uncertainties, the observed limits corresponding to the WIMP-parton cross 
section obtained from the — lctheory lines in figure 4 are shown as thin dotted lines. The latter 
limits are conservative because they also include theoretical uncertainties. The ATLAS limits are 
for the four light flavours assuming equal coupling strengths for all quark flavours to the WIMPs. 
For comparison, 90% CL limits from the SIMPLE [73], Picasso [74], CDF [19], and CMS [21] 
experiments are shown. 



of WIMPs [13, 15]. This is shown in figure 7 where the limits on vector and axial- vector 
interactions are translated into upper limits on the annihilation rate of WIMPs to the four 
light quark flavours. The annihilation rate is denned as the product of cross section a and 
relative velocity v, averaged over the dark matter velocity distribution ((a v)). Equations 
(10) and (11) of ref. [15] are used to calculate the annihilation rates shown in figure 7. For 
comparison, limits on annihilation to bb from Galactic high-energy gamma-ray observations 
by the Fermi-LAT experiment [75] are also shown. The Fermi-LAT values are for Majorana 
fermions and are therefore scaled up by a factor of two for comparison with the ATLAS 
limits for Dirac fermions (see for example the description of equation (34) of ref. [76] for an 
explanation of the factor of two) . Gamma-ray spectra and yields from WIMPs annihilating 
to bb, where photons are produced in the hadronisation of the quarks, are expected to be 
very similar to those from WIMPs annihilating to lighter quarks [77, 78]. In this sense the 
ATLAS and Fermi-LAT limits can be compared to each other. The figure also demonstrates 
the complementarity between the two approaches. The Fermi-LAT experiment is equally 
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Figure 7. Inferred ATLAS 95% CL limits on WIMP annihilation rates (a v) versus mass m x . 
(a v) is calculated as in ref. [15]. The thick solid lines are the observed limits excluding theoretical 
uncertainties. The observed limits corresponding to the WIMP-parton cross section obtained from 
the — lothcory lines in figure 4 are shown as thin dotted lines. The latter limits are conservative 
because they also include theoretical uncertainties. The ATLAS limits are for the four light quark 
flavours assuming equal coupling strengths for all quark flavours to the WIMPs. For comparison, 
high-energy gamma-ray limits from observations of Galactic satellite galaxies with the Fermi-LAT 
experiment [75] for Majorana WIMPs are shown. The Fermi-LAT limits are scaled up by a factor 
of two to make them comparable to the ATLAS Dirac WIMP limits. All limits shown here assume 
100% branching fractions of WIMPs annihilating to quarks. The horizontal dashed line indicates 
the value required for WIMPs to make up the relic abundance set by the WMAP measurement. 



sensitive to annihilation to light and heavy quarks, whereas ATLAS probes mostly WIMP 
couplings to lighter quarks and sets cross-section limits that are superior at WIMP masses 
below 10 GeV for vector couplings and below about 100 GeV for axial-vector couplings. At 
these low WIMP masses, the ATLAS limits are below the value needed for WIMPs to make 
up the cold dark matter abundance (labelled Thermal relic value in figure 7), assuming 
WIMPs have annihilated exclusively via the particular operator to SM quarks while they 
were in thermal equilibrium in the early universe. In this case WIMPs would result in 
relic densities that are too large and hence incompatible with the WMAP measurements. 
For masses of m x > 200 GeV the ATLAS sensitivity worsens substantially compared to the 
Fermi-LAT one. This will improve when the LHC starts operation at higher centre-of-mass 
energies in the future. 
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The value of using an effective field theory approach to WIMP-SM particle coupling is 
that only two parameters, M* and m x , are needed to describe WIMP pair production at 
the LHC, WIMP-nucleon scattering measured by direct-detection experiments, and WIMP 
annihilation measured by indirect-detection experiments. The complementarity between 
the different experimental approaches can hence be explored under a number of important 
assumptions: the effective field theory must be valid, WIMPs must interact with SM quarks 
or gluons exclusively via only one of the operators of the effective field theory (since a mix 
of operators with potential interference effects is not considered here), and the interactions 
must be flavour-universal for the four light quarks. In the future, should there be a WIMP 
signal in at least one of the experiments from these various fields, the effective-operator 
approach would allow important tests of the underlying physics by probing all the available 
experimental data. 

7 Summary 

A search for physics beyond the Standard Model is presented in events with a high-energy 
jet and missing transverse momentum. The search uses the full 2011 pp LHC dataset 
recorded with the ATLAS detector at a centre-of-mass energy of i/s = 7 TeV. The data 
correspond to an integrated luminosity of 4.7 fb . 

Four overlapping signal regions are defined for the search. They require a high-energy 
jet and missing transverse momentum of at least 120, 220, 350 and 500 GeV, with at most 
one additional jet not aligned with the direction of E™ 1SS (to suppress multijet background). 
In all cases the events are required to contain no identified electrons or muons. The 
dominant Standard Model backgrounds from Z and W plus jet production, where the 
boson decays to a final state that includes 1-2 neutrinos, are determined using data control 
regions with correction and transfer factors determined from data and simulations. This 
technique allows precise estimates of the SM contributions to monojet final states, which 
is reflected in a small total uncertainty of 3.2% for the background prediction in the high- 
statistics signal region SRI. 

In each of the four signal regions, agreement is found between the Standard Model 
predictions and the data. Upper limits are set at 95% CL on the visible cross section of 
any non-SM contribution to the signal regions. These limits range from 1.92 pb in the first 
signal region to 7fb in the fourth signal region. To allow comparisons with the results 
of other experiments, 90% CL limits are also provided. The cross-section upper limits 
are interpreted in terms of limits on the model parameters of two BSM physics scenarios. 
For ADD, a model of large extra spatial dimensions, lower limits are set on the (4 + n)- 
dimensional Planck scale Md of 4.17 (2.51) TeV for n = 2 (6) extra dimensions at LO and 
4.37 (2.53) TeV at NLO. In a second scenario an effective field theory is used to derive 
limits on a mass suppression scale M* for pair production of WIMP dark matter particles. 
Within this approach the ATLAS limits can be converted to limits on WIMP-nucleon 
scattering and WIMP annihilation cross sections. Assuming the effective field theory is 
valid, that WIMPs interact with SM quarks or gluons, and that they can be pair-produced 
at the LHC, some of the limits are competitive with or substantially better than limits set 
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by direct and indirect dark matter detection experiments, in particular at small WIMP 
masses of m x < 10 GeV. 
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